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Abstract.  The  smoothing  method  is  used  to  find  the  first  two  moments,  i.e.,  the  mean 
and  the  two-point  two-time  correlation  function,  of  the  field  scattered  by  a  rough  surface. 

The  results  are  expressed  in  terms  of  a  reflection  coefficient  and  a  differential  scattering 
coefficient.  They  are  compared  with  those  found  by  several  other  methods. 

is 

1.  Introduction 

An  acoustic  wave  which  hits  a  rough  surface  produces  a  mean  or  coherent  reflected  field  and  a  fluctuating 
or  partially  coherent  scattered  field.  These  fields  can  be  calculated  by  the  regular  perturbation  or  Born 
expansion  method,  see  t.g.,  [1],  bnt  the  results  diverge  at  grazing  incidence.  However  they  can  be  modified 
or  renorm alised  to  remain  finite  (1)  by  writing  them  in  forms  suggested  by  Twersky’s  self-consistent  field 
method  [2,  3).  We  shall  show  that  these  same  non-divergent  results  can  be  obtained  directly  by  a  different 
perturbation  procedure,  known  as  the  smoothing  method. 

The  smoothing  method  has  been  used  to  find  the  coherent  field  for  certain  kinds  of  rough  surfaces  by 
Wenzel  [4],  DeSanto,  and  others.  (See  DeSanto’s  review  [5].)  Here  we  apply  it  to-  more  kinds  of  rough 
surfaces,  which  can  be  either  moving  or  at  rest,  and  we  calculate  both  the  coherent  field  and  the  two-point 
two-time  correlation  function  of  the  field.  Then  we  show  that  the  results  are  the  same  as  the  renormal- 
ixed  ones  in  [1),  which  were  compared  there  with  Twersky’s  [2,  3]  results  for  embossed  surfaces.  Thus  the 
present  calculation  elucidates  the  relationship  among  the  results  of  the  regular  perturbation  or  Born  expan¬ 
sion  method,  the  renormalized  regular  perturbation  method,  Twersky’s  self-consistent  field  method  and  the 
smoothing  method. 

3.  The  Integral  Equation  for  the  Scattered  Field 

Suppose  that  a  wave  with  the  acoustic  velocity  potential  i>(z,  y , z, t)  is  incident  upon  a  rough  surface  S,  and 
that  it  produces  a  scattered  wave  with  potential  <p{z,  y,  z,  t).  Both  potentials  must  satisfy  the  wave  equation 
with  constant  sound  speed  c  above  the  surface  S,  and  f>  must  be  outgoing.  In  addition  +  ifi  must  satisfy 
on  S  a  boundary  condition  which  depends  upon  the  nature  of  the  surface.  We  shall  consider  the  following 
four  boundary  conditions: 

5«(l4  +  v)  =  0  on  z  =  eh  (la') 


(d,  -  ch,d,  -  +  <p)  =  cht  on  z  =  eh 


(15') 


[chd,  +  ca,)(V>  +  v>)  =  0  on  z  =  0  (Ic) 

I 

(3f  +  tehd,)(^>  +  p)  =  0  on  z  =  0  (Id) 

In  all  cases  £  is  a  small  parameter  and  h[x,  y,  t)  is  a  given  function.  Case  a  represents  a  soft  surface,  z  =  eh, 
on  which  the  pressure  vanishes;  ease  4  represents  a  hard  surface  z  =  eh  on  which  the  acoustic  normal  velocity 
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equals  the  surface  normal  velocity,  case  c  represents  a  Bat  surface  z  =  0  with  the  admittance  eh,  and  case  4 
represents  a  Bat  surface  z  =  0  with  the  impedance  eh.  In  case  i  we  shall  not  consider  the  acoustic  radiation 
produced  by  the  motion  of  the  surface,  so  we  shall  omit  the  term  eh,  on  the  right  side  of  (lb'). 

Before  proceeding  we  shall  simplify  (la')  and  the  homogeneous  form  of  (lb')  by  expanding  +  V  in  a 
Taylor  series  about  z  =  0  to  obtain 

,lcl 

[1  +  thd,  +  — — -d]  +  0(es)]3,(tb  +  p)  =  0  on  z  =  0  (la) 

* 

[<?,  +  e(h<?J  -  htdx  -  hyd,)  +  e5(— dj  -  hhtdtd,  -  hh,<?#<9,)  +  0(e3)](p  +  i/>)  =  0  on  z  =  0  (lb) 

Now  all  four  boundary  conditions  (la)-(ld)  apply  on  the  plane  z  =  0.  This  makes  it  convenient  to  introduce 
the  Fourier  transform  pairs: 

F(»)  =  (2*)-*J  /(p)e->“dp  (2) 

/(p)  =  j  (3) 

Here  p  =  (x,y.t)  and  s  =  (a,/3,w)  are  three  component  vectors  in  physical  space-time  and  in  wavenumber- 
frequency  space,  respectively.  We  shall  always  use  the  corresponding  capital  letter  to  denote  the  transform 
of  a  function  denoted  by  a  small  letter. 

We  now  apply  the  transform  (2)  to  the  wave  equation  for  p  and  find  that  $(*,  *)  satisfies  the  ordinary 
differential  equation  +  fc5(s)<b  =  0  where  fc5(a)  =  uiJ/cJ  -  oJ  —  /3J.  Since  ♦  must  be  outgoing  at  *  =  +oo 
when  k  is  real,  and  must  vanish  there  when  k  is  imaginary,  4  must  be  of  the  form 

♦(*,*)  =  d.(s)e-“<‘)'  (4) 

Here  k(s)  is  defined  by 

fc(s)  =  («/«:)(l-(a,+/J*)(<-/c)-*)»  ,  +  I J*  (5) 

*(«)  *  -»'(o*  +  j9*-(w/e)*]i  ,  (w/e)’  <a*  +  0* 

The  scattered  amplitude  d,(s)  is  unknown,  and  is  to  be  found.  To  find  it  we  first  write  t(s,  z),  the  transform 
of  the  incident  field,  as 

•  (.,*)  =  d<(s)e*<‘>*  (6) 

The  incident  amplitude  i4,(s)  is  assumed  to  be  knowu.  We  now  apply  the  transform  (2)  to  each  of  the 
boundary  conditions  (la)  -  (Id)  and  use  (4)  and  (6)  in  the  resulting  equations.  Is  each  case  we  obtain  an 
integral  equation  for  d,(a)  of  the  form 

K[-t)A.  =  ±K[e)Ai  (7) 

The  linear  operator  K{e)  m  (7),  which  depends  upon  c,  can  be  written  as 

K(e)  =  I  +  cKi  +  eiKt  +0(t*)  (8) 

Here  I  is  the  identity  while  Ky  and  Kt  are  integral  operators  of  the  forms 

KyA  =  /  H(.  -  s')m(s,  s')d(s')d#'  (9) 


**A  =  -i-JJ  B(a  -  .')»(«'  -  .'')[^]fcV)d(.'')d."ds'  (10o) 


•  pp  "C.  j  !  V  . 

U!  otrit-ut.  tr.  :  .  : . 

HATTHS#  J.  KE1. 

Chlsf,  Technlcil  1  nf  onnat  1  on  PI  vis  t  on 


=  \(j  {**(«")  -2[K  /a)*  -a'Q"-ff0"]}A(*")dt"'U'  (104) 

We  recall  that  H  ia  the  transform  of  h,  while  m  and  the  choice  of  sign  in  (7)  are  given  in  Table  I.  For  the 
admittance  and  impedance  boundary  conditions,  cases  e  and  d,  the  operator  A'j  and  the  0(c3)  term  in  (8) 
are  both  tero. 


- (SI - (HI 

(') 

131 

Parameter 

Soft  Hard 

Admittance 

Impedance 

± 

+ 

+ 

- 

m(s,s') 

iu'k(»')/u  sfc-»(a)  [(*)*  -  oa'  -  /J/}'] 

u i' /ck(s) 

w/eJt(a') 

Table  I.  The  sign  and  the  function  m(s,s'). 


3.  Application  of  the  Smoothing  Method 

Let  h  and  A,  be  random  functions  statistically  independent  of  one  smother.  Then  the  solution  A.  of  the 
integral  equation  (7)  is  random  also.  To  calculate  its  first  two  moments  we  shall  use  the  smoothing  method, 
first  used  to  treat  waves  in  random  media  and  similar  problems  by  Bonrret  [6]  and  Keller  [7].  It  is  convenient 
to  write  (/)  to  denote  the  mean  value  with  respect  to  h  of  any  random  function  /  that  depends  upon  h, 
and  /'  =  /-(/)  to  denote  the  corresponding  fluctuating  part  of  /.  Then  {/')  =  0.  If  /  depends  upon  both 
h  and  Ai,  then  {/)  still  depends  upon  A,-,  so  it  is  still  random.  We  denote  the  average  of  /  with  respect  to 


both  h  and  Ai  by  ({/)). 

We  begin  by  using  (8)  in  (7)  to  get 

(I  -  tKt.  +  e’.KsjA,  =  ±(/  +  eK,  +  t*Kt)Ai  +  0(e»)  (11) 

Now  we  take  the  mean  value  of  (11)  with  respect  to  A  to  obtain 

(/  -  e(Ki)  +  t'(Ki)){A.)  -  e(K(A't)  +  e'(K'tA',)  =  ±(I  +  e(Kt)  +  e*(Kt))Ai  +  0(e»)  (12) 

Next  we  subtract  (12)  from  (11),  keeping  terms  through  0(e),  to  get 

A'.  -  t(K[(A.)  +  (ffj)A',  -  K[A'.  +  <K;a'.»  =  ±cK(Ai  +  0(t *)  (13) 

Upon  solving  (13)  for  A',  up  to  0(e),  we  find 

A'.=eK[(A.)±eK'lAi  +  0(ei)  (14) 

Finally  we  use  (14)  for  A',  in  (12)  and  obtain 

(/  -  t(Ki)  -  eV/fitfi)  +  *’(A-.>](A.)  =  ±[1  +  e(JCt)  +  t'(K[K()  +  e*(JT,)]A(  +  0(e»)  (15) 


This  is  the  desired  smoothed  or  non-random  equation  for  (A,),  and  (14)  is  an  expression  for  A'  in  terms  at 
(A,).  The  operators  (Kt)  and  (Kt)  in  (15)  are  given  respectively  by  (9)  with  H  replaced  by  (H)  and  by 
(10)  with  HH  replaced  by  (HH).  The  operator  (K(K[)  is  given  by 

(K(K()A  =  J (H'(,  -  s')£T'(s  -  ."))m(,,s')m(s\s")A(s")ds'ds"  (Id) 
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The  integral  equation  (15)  simplifies  considerably  when  the  function  h  is  statistically  second  order 
stationary.  Then  the  mean  of  h  is  a  constant  h„,  i.e.,  (h)  =  ha,  and  the  two-point  two-time  correlation 
function  of  h  depends  only  upon  the  difference  between  its  arguments,  i.e.  (h(p  +  p')h(p'))  =  r(p)  +  hj. 
Consequently  the  first  two  moments  of  the  transform  H  axe  given  by  (ff(s))  =  h„<5(a)  and  {/f(s)//(s'))  = 
R(a)f(a  +  a')  +  h’i(s)6(j').  Here  R{a),  the  spectral  power  density  of  h(p).  is  the  Fourier  transform  of  the 
auto-corTelation  function  r(p),  so  it  is  a  positive,  even,  real- valued  function  of  a.  When  the  expressions  for 
(H)  and  (H(a)H(a'))  are  used  in  the  average  operators,  they  simplify  to, 

{Ki)A  =  h0m(a,a)A(a)  (17) 

{K{K\)A  =  [J  R(a  -  s')m(s,s')m(s',s)ds']d(s)  (18) 

(Ki)A  =  -\k\a)  [V  +  h*]A(s)  (19a, fc) 

Here  V  =  /  R(a)da  is  the  variance  of  h.  Thus  for  stationary  processes,  these  three  operators  are  multi¬ 
plicative,  and  the  smoothed  equation  (15)  is  just  a  linear  algebraic  equation  for  {A.).  In  cases  c  and  d, 


(*.)  =  0 

4.  Solution  for  the  Mean  Amplitude 

In  general  the  solution  of  (15)  for  {A,)  can  be  written 

(A.)  =  MAt  +  0(e3)  (20) 

The  linear  operator  M,  which  occurs  in  (20),  is  not  random,  and  most  be  found  by  solving  the  integral 
equation  (15).  However,  when  h(p)  is  second  order  stationary,  the  solution  of  (15)  is 

(A.(s))  =  C(a)A<(.)-i-0(e‘)  (21) 

Here  C(s)  is  just  a  scalar  function  which  we  call  the  reflection  coefficient.  It  is  given  by  . 

C(a)  =  ±[1  +  e(JCi)  +  e*(K[K[)  +  c’(ff ,)]/[!  -  t{Kt)  -  t'{K[K[)  +  e*{Kt)\  +  0(e5)  (22) 

=  ±(1  +  <?(*)]/(l  -  Q(s))  +  0(«*) 

The  sign  is  shown  in  Table  I  and  <J(»)  is  defined  by 

<?(*)  =  ^eh.m(s,  s)  +  e’  J  fl(a  -  «')m{s,  a')m(s,  »')da'j  |l  -  yi’(a)  (V  +  hj)|  (23a, 6) 

Q(a)  =  eh0m{»,a)+ei  J  R{*  -  »')m(j,  «')<*«'  (23c,  d) 

When  (21)  is  averaged  with  respect  to  Aj,  it  yields 

«A.(a)))  =  <7(.)<A,) +  <?(«’)  (24) 


Here  the  average  of  A,  pertains  to  both  h  and  A<,  while  (A<)  means  the  average  of  A,,  which  is  independent 
of  h. 

The  second  form  of  C(s)  in  (22)  is  obtained  from  the  first  form  by  dividing  numerator  and  denominator 
by  1  +  c*(K j),  following  Twersky  [8).  If  the  denominator  in  (23a, b)  is  replaced  by  unity,  which  corresponds 
to  an  O(c’)  change  in  C(«),  the  resulting  C(«)  is  identical  with  equation  (42)  of  [1],  which  is  the  renormalised 
Bom  reflection  coefficient. 
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6.  The  Second  Moment  of  the  Field 

Let  us  consider  the  total  acoustic  potential  i/>  +  <p  » t  two  points  with  separation  x  =  (x,  y,  z,  t)  and  midpoint 
X  =  (X,  V,  Z,T).  The  second  moment  or  two-point  two-time  correlation  function  of  the  acoustic  potential 
at  these  points  is  -j(x,  X)  =  (([*  +  *p\(X  +  +  v\(X  —  f ))).  We  denote  its  Fourier  transform  r(a,q,z,  Z) 

where  s  =  (a,/},ui)  corresponds  to  (x,y,t)  and  q  =  (a,  6,  w)  corresponds  to  (X,  Y,  T),  and  we  write  it  as 

T[a,q,z,Z)  =  «[*  +  *](i  +s,Z  +  1)(*  +  *](|  -  s,Z  -  f  ))>  (25) 

By  using  (4)  and  (6)  for  *  and  *  we  can  express  T  in  terms  of  the  four  correlations  Bu(a,q)  = 

UMt  +  «Mi(f  ~  •)».  =  ((AtA.)),  Bti  =  {(A.Ai))  and  £„(.,?)  =  ((A.(*  +  >)A.(\ -  •))).  See  [I]. 

Each  of  the  four  correlations  above  can  be  evaluated  by  writing  A,  =  (X,)  +  A\  with  (A,)  given  by  (20) 
and  A',  given  by  (14).  In  particular  the  auto-correlation  of  the  scattered  amplitude  4,  becomes 

*..(•.«)  =  «*.(§  +  •)><*.(}  ~  •)»  +  ~  *))>  +  0(e5)  (26) 


=  {[MA^l  +  »)1M^,)(|  -  •))  +  **({IJrjr  V  ±  I)A,}C-  +  s)[Xj(Af  ±  /)A.](|  -  s)»  +  0(es) 
When  h  is  second  order  stationary,  M  becomes  multiplication  by  C(«),  and  (26)  simplifies  to 


B..(>,q)  =  C(|  +  s)C(|  -s)Bii[»,q)  +  4c1  j 


J  f  R(t-  s')rn(f  +  s,  *  +  s')m(  J  -  s,  f  -  a')Bn(a',q) 

J  [1 -0(f +  *')][! -Q(f-e')l  +°le) 


The  auto-corTelation  B„  simplifies  even  more  when  the  incident  field  *(x,  y,  z,t)  is  statistically  second 
order  stationary  in  x,y  and  t.  Then  Bu  is  just 

Bu{a,q)  =  /<(«)%)  (28) 

where  I((a)  is  the  intensity  of  the  incident  field.  With  (28)  used  in  it,  (27)  becomes 

B..(s,q)  =  T.{a)S(q)  +  0(e*)  (29) 

The  function  /,(*)  in  (29),  which  is  the  intensity  of  the  outgoing  field  in  the  direction  s,  is  defined  by 

/.(.)  =  |C(s)|>A(.)  +  (fc-*(s)|  j  <r(e,*')A(*')da'  +  0(«»)  (30) 

The  coefficient  <r(«,  s')  in  (30)  can  be  identified  as  the  differential  scattering  coefficient  of  the  surface.  It  is 
defined  by 

(  J  -  |1  -  Q(,')|»  (31) 

With  the  denominator  in  (23a, b)  replaced  by  unity,  this  is  exactly  the  result  for  <r  given  in  equation  (44)  of 
[1],  and  called  there  the  renormalized  Bora  approximation. 


0.  Conclusion 

By  using  the  smoothing  method,  we  have  calculated  the  first  two  momenta  of  the  field  produced  when  a 
possibly  random  incident  field  hits  a  random  slightly  rough  surface.  The  first  moment  of  the  field,  called  the 
coherent  field,  consists  of  an  incident  and  a  reflected  field.  The  reflected  field  is  determined  by  a  reflection 
coefficient  C(«),  which  depends  upon  the  first  two  moments  of  the  surface  roughness  h.  The  second  moment 
of  the  field  involves  a  differential  scattering  coefficient  <r(», »')  which  also  depends  upon  the  second  moment 
of  K 


1 


Our  results  for  C(s)  and  <7(3,  a')  are  the  same  as  those  wc  obtained  before  [1]  by  renormalixing  the 
divergent  results  given  by  the  Born  approximation.  That  renormalization  was  achieved  by  writing  C  and  o 
in  the  forms  obtained  by  Twersky  [2,3J.  Wc  also  showed  in  [1]  that  for  an  embossed  plane,  our  results  are  the 
same  as  his  except  for  one  difference.  We  had  the  second  Bom  approximation  to  the  differential  scattering 
amplitude  of  a  single  boss  instead  of  the  exact  amplitude  which  occurs  in  his  theory.  Tliis  same  comparison 
applys  to  the  results  of  the  smoothing  method,  as  used  here. 

Wc  conclude  that  three  different  methods  yield  the  same  results  for  the  C  and  a  of  an  embossed 
plane,  provided  that  each  boss  is  small  enough,  with  a  small  enough  slope,  so  that  its  differential  scattering 
amplitude  is  well  approximated  by  its  second  Bom  approximation.  The  three  methods  are  Twcrsky's  method, 
the  renormalized  Born  approximation  and  the  smoothing  method  as  used  here.  When  the  bosses  or  their 
slopes  are  not  small  enough,  Twersky's  results  are  better  than  the  other  two.  On  the  other  hand,  if  the 
surface  is  not  an  embossed  plane,  Twersky's  method  does  not  apply  but  the  other  two  methods  do. 
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